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O ■ Abstract 

Pseudo-scalar and vector meson correlation functions were calculated 
' at temperatures below and above the deconfinement transition using 

rjj^ O(a) improved Wilson fermions in quenched lattice QCD. The spectral 

functions were reconstructed from the correlator given only at discrete 
| points in Euclidean time by means of the Maximum Entropy Method 

without a priori assumptions on the spectral shape. The deviation of 
the spectral functions from the one of a freely propagating quark anti- 
quark pair is examined above the critical temperature. At 1.5 and 3 T< 
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C*~> • the vector spectral function yields an enhancement of the dilepton rate 

over the Born rate in the energy interval 4 < uj/T < 8 and a sharp drop 
O ' at energies below 2-3 T. 
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1 Introduction 
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The thermal dilepton spectrum is accessible in heavy ion collisions and pro- 
vides a key observable to study thermal properties of the medium at high 
density and temperature |], ||. Thermally induced changes of the dilepton 
spectrum at low energies |§ are expected to be influenced by non-perturbative 
in-medium modifications of the quark anti-quark interactions. Furthermore 
such effects are directly related to changes of the spectral shape in the vector 
meson channel. In two-flavor QCD the differential dilepton rate is connected 
to vector spectral function cry (to, p, T) through the relation 

dN a = dW 5a 2 1 ( - T] m 

d*xd*p ~ dud^p 27tt 2 uj 2 {e"/ T - 1) 1} ' {1) 

At the same time the spectral densities are related to the Euclidean cor- 
relation functions of hadronic currents Jh{t,x) = ^>(t, x) Th ip{r,x) with 
Th G [1, 75, 7 M , . . .] through the integral equation 
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Although hadronic correlators can be calculated numerically in the framework 
of lattice regularized QCD for discrete Euclidean times r e [0, 1/T], the in- 
version of the integral equation still remains an ill-posed problem. This situa- 
tion improved considerably after the implementation of the Maximum Entropy 
Method (MEM) for lattice QCD [|j], which allows to obtain the most probable 
spectral function without a priori assumptions on the spectral shape. 

Most of the previous applications of MEM concentrated on the analysis of 
hadronic properties like masses and decay width at zero temperature, where the 
pole masses of the ground f|, |] and even excited states || could be obtained 
with quite satisfactory precision. The great challenge, however, is to apply the 
method at non-zero temperature, where very little about the spectral shape is 
known so far. 



2 Thermal Meson Spectral Functions 

At finite temperature excited states and continuum- like contributions are ex- 
pected to gain in influence on the spectral shape. Thus it is mandatory to 
verify that the method allows to reconstruct such spectral functions correctly. 
A first step in this direction was the reconstruction of the spectral shape for free 
massless quarks in the (pseudo-)scalar channel Op^ 6 = 3/8 n~ 2 u) 2 tanh(u;/4T), 
known from leading order perturbation theory. In this case the relation (0) was 
used to obtain discrete values of the correlator in Euclidean time j5|. On the 
other hand such a correlator is directly calculable on the lattice 0, which is 
shown in figure 1(a) for three different lattice sizes. At small time separations 
the finite lattice cut-off effects are visible compared to the free continuum cor- 
relation function indicated as solid line. Nevertheless, the spectral shape could 
be reconstructed almost perfectly already for temporal extents N T = 8 . . .16 
(see fig. 1(b)) by introducing a lattice version of the integration kernel in equa- 
tion (Q) for vanishing momentum 

G h (t,T)= dLua v (u,T)K L (T,u,N T ) , (3) 
Jo 

where Kl(t, uj, N t ) is the finite lattice approximation of the continuum kernel 

2uj n ^ exp(-i27mrT) 

K l (t,uj,N t ) = — ^ — . 4 

T ~ (2N T sm(mT/N T )y + (lu/T) 2 

Since the applicability of MEM at finite temperature has been verified, one 
can proceed to analyze meson spectral functions below and above the criti- 
cal temperature of the deconfinement transition. We have calculated pseudo- 
scalar and vector meson correlation functions on large isotropic lattices of sizes 




Figure 1: Free scalar correlators (a) and reconstructed spectral shape (b) 



(24 — 64) 3 x 16 in a temperature range 0.4 to 3.0 T c using 0(a) improved Wil- 
son fermions || in quenched lattice QCD. Below T c up to five quark mass 
values were chosen in order to obtain a reliable extrapolation to the chiral 
limit (m q — > 0). Above the critical temperature the simulations could be per- 
formed directly at approximately zero quark mass, since there are no longer 
massless Goldstone modes present ||. 

Compared to the sharp (^-function like peaks observed in the spectral func- 
tions at T = a broadening of the peaks corresponding to particle poles 
and a reduction in height can be observed at finite temperature. Figure 2 
shows an example at 0.4 T c , where the error bars indicate the uncertainty in 
the given energy region calculated from the covariance matrix of the spectral 
function |IJ. The ground state 
masses extracted from the spectral 
functions at 0.4 T c agree very well 
with zero temperature data ||11||. In- 
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creasing the temperature towards T c 
the vector meson spectral function 
is investigated at fixed pion mass. 
The ground state peak experiences 
a further broadening and is slightly 
shifted towards higher energies. It 
remains to be clarified to what ex- 
tent this effect is induced by the re- 
stricted temporal dimension and lim- 
ited statistics. 

Above the critical temperature only a gradual approach towards the free 
correlation function is observed in the (pseudo-) scalar channel, while the vector 



10 



20 



30 



40 



50 



60 



Figure 2: Meson spectral functions for 
a ratio of m^/rrip = 0.6 at 0.4 T c 
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Figure 3: Meson correlators (a) and reconstructed spectral functions (b) 



meson correlator is much closer to the free quark behavior already at 1.5 T c 
(see fig. 3(a)). Analyzing directly the ratio of the correlation function in the 
vector channel to the corresponding free curve calculated for the same lattice 
size yields a deviation of about 10 %. The free correlator is approached from 
above, while simple quasi-particle pictures lead to a reduction of Gy(r,T) 
relative to G f v ree (r,T) Q, pf. The deviat ions from the free quark behavior are 
also evident from the reconstructed spectral functions in figure 3(b). While the 
second enhancement over the free curve visible in fig. 2 and 3(b) at energies 
larger than uj/T = 16 is presumably a lattice artifact due to the heavy fermion 
doublers of the Wilson fermion formulation ||, the peak at about uo/T = 6 
shows the actual physical effect of persisting interactions between quarks and 
gluons in both channels up to 3 T c . Compared to the zero temperature results 
[|, H, |6| the peaks are rather broad and less pronounced in height, but might 
sharpen by increasing the statistics and the lattice size @]. 

3 Thermal Dilepton Rates 

Once the most probable vector meson spectral function is reconstructed from 
the correlation function in Euclidean times by means of the Maximum Entropy 
Method, the remaining step to obtain the differential dilepton rate is straight 
forward. The spectral functions shown in figure 4(a) were obtained from corre- 
lators on the largest lattice (64 3 x 16) with approximately zero quark mass at 
temperatures 1.5 and 3.0 T c ||. In addition to the purely statistical error-band 
obtained from a Jackknife analysis, the uncertainty of the uniqueness of the 
result incorporated in MEM is shown in the insertion as error bars for four 
energy intervals. This leads to the conclusion that the data sample is large 
enough to yield statistically significant results. 




Figure 4: Vector spectral function (a) and differential dilepton rate (b) 



The differential dilepton rate can now be calculated in employing relation 
([I]) projected onto zero momentum. Figure 4(b) illustrates the correspond- 
ing dilepton rates for both temperatures. In the high temperature limit the 
leading order perturbative result can be obtained by using free massless quark 
propagators in the calculation of the current-current correlation functions. For 
vanishing momentum the spectral function in the vector channel is given by 
atj ee = 3/4 7t~ 2 cg> 2 tanh(o;/4T) , which leads to the reduction of the differential 
dilepton rate to the Born rate 

djy Born _ ha 2 1 

du;d 3 p ^ = 0) = 36vr 4 (e-/ 2T + l) 2 ' (5) 

The broad peak visible in the spectral functions at around uj/T = 5 — 6 results 
in an enhancement of the respective dilepton rates in the region uj/T = 4 — 8 
compared to the Born rate, which seems to scale with the temperature. For 
even larger energies the close agreement of the differential dilepton spectrum 
with the Born rate is evident. 

The most striking feature of the observed dilepton rate is certainly the 
sharp drop at energies below uo/T = 2 — 3. This effect is in contrast to hard 
thermal loop resummed perturbation theory |T3j as well as to 2-loop perturba- 



tive calculations which lead to a divergent vector spectral function in the 
limit uj — > 0. If this suppressions persists in future investigations closer to the 
critical temperature, there would be no thermal contribution the dilepton rate 
at small energies during the expansion of the hot medium created in heavy ion 
collisions. This is consistent with the present results at SPS energies |15[] and 
might be observable in the dilepton spectra at RHIC or LHC energies. 
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